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ABSTRACT 

We present ULTRACAM photometry of ES Get, an ultracompact binary with a 620s 
orbital period. The mass transfer in systems such as this one is thought to be driven 
by gravitational radiation, which causes the binary to evolve to longer periods since 
the semi-degenerate donor star expands in size as it loses mass. We supplement these 
ULTRACAM+WHT data with observations made with smaller telescopes around the 
world over a nine year baseline. All of the observations show variation on the orbital 
period, and by timing this variation we track the period evolution of this system. We do 
not detect any significant departure from a linear ephemeris, implying a donor star that 
is of small mass and close to a fully degenerate state. This finding favours the double 
white dwarf formation channel for this AM CVn star. An alternative explanation is 
that the system is in the relatively short-lived phase in which the mass transfer rate 
climbs towards its long-term value. 
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1 INTRODUCTION 

Ultracompact binaries with periods of the order of tens 
of minutes or less have attracted much attention in re- 
cent years. These systems consist of a primary white dwarf 
with a companion star that is also at least partially de- 
generate. Close double-degenerate binaries are one of the 

J iroposed progenitor populations of Type la superno vae 
Kotakll2008l : lDi Stefandl201ol : ICilfanov fc Bogdanll20ld ). as 
well as the rec ently proposed sub- luminous '.la' supernovae 
jBildsten et al . 2007; Kasliwal ct al. 2010.). They are of in- 
terest from a binary formation and evolution point of view, 
with the short periods implying at least one common en- 
velope phase in the history of the binary, and the chemical 
composition suggesting helium white dwarfs, helium stars 
or cataclysmic varia bles (CVs) with evolved sec ondaries as 
possible progenitors (jNelemans et al.ll200ll . [2OI0I ). The mass 
transfer in these systems is thought to be driven by angular 
momentum loss as a result of gravitational radiation. These 
sources are predicted to be a mong the strongest gr avita- 
tional wave sources in the sky (|Nelemans et al.ll2004l ). and 



are the only class of binary with examples already known 
which will be detectable by the gravitational wave obse rva- 
tory LISA l|Stroeer fc Vecchidl2006l : iRoelofs et al.ll2007l ). 

The two sy stems with the shor test known periods are 
HM Cnc (324s, iRoelofs et al.ll201Gl ) and V407 Vul (569s, 
iHaberl fc Motchlll995l ')." The exact nature of these systems 
is unclear: the t wo leading mod els are the unipolar induc- 
tion (UI) mod el (|Wu et al. 2002) and the direct-impact ac - 
cretion model jMarsh fc Steeghs.2002. : iRoelofs et al.ll201Gl ). 
There are also 23 known systems with longer periods rang- 
ing from 620 to 3906s (the AM Canum Venaticorum stars, 
AM CVns; see Solhoim 2010 for a recent review). The na- 
tures of these systems are better established: the spectra 
show an absence of hydrogen and the presence of helium 
lines, many of which are the double-peaked e mission lines 
characteristic of sources accreting via a disc l|Marsh|[l999l : 
[Moralcs-Rueda et al. 2003). These systems are clearly ac- 
creting, and are the helium equivalent of the cataclysmic 
variable (CV) stars. 

Gravitational radiation has a huge influence on these 
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systems, driving the evolution and determining the or- 
bital period distribution, luminosities and numbers. Mea- 
surements of the time derivative of the orbital period P 
for V407 Vul and HM Cnc ha ve shown the peri o ds in 
these systems to b e decreasing i Hakala et alj |2003| . 12004 
I Ramsay et al ] I2OO5I : IStrohmaveill2005l l. This is consistent 
with the UI model, and contrary to what might be ex- 
pected for accretion. In an accreting double-degenerate sys- 
tem, the gravitational radiation should drive the binary to- 
wards longer orbital periods. However, if the mass trans- 
fer in an accreting double-degenerate system is significantly 
lower than its equilibrium value, due to either some non- 
secular process l|Marsh fc Nelemansll2005|') or the binary be- 
ing in its mass transfer turn - on phase~(IWillems fc Kalogeral 
I2OO5I : iD'Antona erahl l2006l : iDelove fc Taamll2006D . then a 
decreasing binary period is possible. The counter-argument 
to this is that such phases are expected to be short, al- 
though lOelove fc TaamI l|2006l) calculated that the early con- 
tact phase can last much longer (10"^ - 10® yr) than originally 
thought. 

To date, no period derivative has been measured in any 
of the longer period AM CVn stars. Since these are unam- 
biguously accreting binaries, then we would expect to detect 
an increasing period in these objects. Additionally, it is not 
certain that gravitational radiation is the dominant angu- 
lar momentum loss mechanism for the ultra-compact binary 
stars: some other mechanism, such as the magnetic brak- 
ing observed in cataclysmic variables, could be operating to 
drive th e evolution at a higher r ate than assumed in current 
models ({Farmer fc Roelofa 2010f ). These astrophysical issues 
will in the future be important for the use of these systems 
for the verification of LISA a s they can make the diffe rence 
between detectability or not (jStroeer fc Vecchioll2006f ). 

In order to address this issue we have conducted a 
long-term timing study of the AM CVn star ES Ceti 
l| Warner fc Woudj [iooa ) . This system shows a 620s opti- 
cal modulation which is c omplex in structure and varies 
on a night-by-night basis ()Espaillat et al.ll2005l ). Spectro- 
scopic observations have confirmed that this is the orbital 
period, and also show the double-peaked helium emission 
lines which imply the presence of an accretion disc in this 
system (Steeghs, in prep). ES Get is the ideal subject for 
a timing study, since it is the shortest period ultracompact 
binary population after V407 Vul and HM Cnc. It poten- 
tially connects these two systems to the rest of the AM 
CVn population: the period of ES Cet is only 51s greater 
than V407 Vul, so ES Cet may also be in the mass transfer 
turn-on phase. If however ES Cet has already evolved on 
to the long-term and stable AM CVn path of lengthening 
period, then we would expect the period evolution to be the 
most rapid of all the systems, since gravitational radiation 
is strongest when masses are large and periods short. 



made using the SDSS it', g' and r' or i' filters. The CCDs 
were windowed in order to obtain exposure times of 2 - 4s, 
depending on conditions. The dead time between exposures 
for ULTRACAM is ~25ms. In 2005 we observed ES Cet five 
times over the period 27 Aug - 01 Sep 2005. Conditions 
were good, with seeing < 1" and good transparency. We ob- 
served ES Cet once more on the night of 07 Jan 2010. These 
observations were affected by poorer seeing and some cloud. 

All of these data were reduced with aperture pho- 
tometry using the ULTRACAM pipeline software. We fol- 
low the same steps for reduction and flux calibration as 
detailed in our previous ULTRACAM publications (e.g. 
ICopperwheat et aLllioiOal lbh. 



2.2 WHT/ACAM 

We made a further observation of ES Cet in excellent condi- 
tions on the night of 22 August 2010 with the WHT us ing the 
perm anently- mounted imaging camera ACAM (Ben n et al.l 
l2008l ). We used a 5 second exposure time and windowed the 
chip to reduce the readout time between exposures to ~0.8s. 
These data were reduced in the same way as the ULTRA- 
CAM data. 



2.3 SAAO and CBA observations 

We supplemented our WHT observations with additional 
data collected with smaller telescopes over a longer time 
period. We observed E S Cet with the Un iversity of Cape 
Town CCD photometer (|0'Donoghuell995l ). mounted on the 
1.9m telescope at the South African Astronomical Observa- 
tory (SAAO). This instrument was used in frame transfer 
mode and the observations were made in white light in or- 
der to maximize the count rate. Typical exposure times vary 
from 10 - 20s. These data were reduced using the DoPhot 



program (|Schechter et al.lll993l : [0'Donoghuelll995l ). We dis- 
carded any data in which conditions, the exposure time or 
the length of the observation means that the photometric 
period is not readily apparent. This left us with 38 nights of 
observation obtained between 18 Oct 2001 and 9 Oct 2009. 
The length of each observation varied, but was typically be- 
tween 5 and 10 orbital cycles. 

We obtained additional nights of observation with the 
Center for Backyard Astrophysics (CBA) telescope network, 
a global network of small telescopes dedicated to photom- 
etry of cataclysmic variable^. Typical exposure times for 
these observations were ~20s. As with the SAAO data, we 
discarded any observations in which the photometric period 
is not obvious. This leaves us with 21 observations obtained 
between 18 Jan 2002 and 1 Jan 2010. 



2 OBSERVATIONS 
2.1 WHT/ULTRACAM 

We observed ES Cet with the h igh speed CCD camera 
ULTRACAM (|Dhillon et al.l 120071 ) mounted on the 4.2m 
William Herschel Telescope (WHT) in 2005 and 2010. UL- 
TRACAM is a triple beam camera and all observations were 



3 ULTRACAM LIGHT CURVES 

We measured the flux from ES Cet in 2005 using our UL- 
TRACAM data, and find the mean magnitudes to be 15.94, 
16.28, 16.63 and 17.23 in u', g' , r' and i' respectively. We 
found the object to be slightly fainter in 2010, with mean 



^ More details concerning the CBA can be found at 
,http: / /cba. phys.coIumbia.edu 
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Figure 1. Phase folded ULTRACAM and ACAM light curves. Due to the significant variations in the light curves from night- to-night, 
we phase fold the data on a night by night basis. The ti'-band light curves are plotted in the left column, the g'-band light curves in the 
middle, and the right column contains either the j'-band (nights 1 and 2) or r'-band (nights 3 to 6) light curves. We omit the u'-band 
light curve for night 6, which was particularly affected by poor weather conditions. A flux offset is applied to each light curve for clarity. 
For the g'-band data we overplot the fits as detailed in Section |4l consisting of three sinusoids at the fundamental orbital frequency and 
its first and second harmonics. 



magnitudes of 16.19, 16.35 and 16.73 in u' , g' and r' . The 
amplitude of the light curve varies between bands, but is 
0.11 - 0.18 on average. 

Using the ephemeris given in Section U we phase- folded 
our ULTRACAM and ACAM data and plot the results in 
Figu re [1] We do no t com bine data from multiple nights, 
since lEspaillat et"al] (|2005l ') showed evidence for some vari- 
ation in the shape of the light curve on different nights. 
We also observe these changes in our data: while each light 



curve is roughly sinusoidal with the minimum light consis- 
tently at the phase of zero, on 28 and 29 August 2005 we 
see a secondary minimum at a phase of ~0.35, which is not 
present in the other three hght curves obtained in 2005. The 
two 2010 hght curves may be intermediate between the two 
states observed in 2005. The timescale over which the light 
curve shape can change is very short: an examination of 
the unfolded hght curves shows significant changes over the 
course of a few orbital cycles. The cause of this variation is 
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Table 1. Our minimum-light timings for ES Get, using the sine 
function fit described in Section|4] We determine the uncertainties 
by fitting a large number of data sets, resampled from the original 
data using the bootstrap method llEfronll979l : lEfron fc Tibshiranil 
For each timing we indicate if it was obtained by SAAO 
(S), CBA (C), ULTRACAM (U) or ACAM (A). 





Cycle 


BMJD 




Cycle 


BMJD 


s 





52200 


98063(3) 


C 


57688 


52615 


08678(6) 


s 


f34 


52201 


94270(3) 


c 


60619 


52636 


12651(3) 


s 


274 


52202 


94762(3) 


S 


60994 


52638 


81849(10) 


s 


433 


52204 


08908(3) 


s 


83192 


52798 


16412(4) 


c 


f2698 


52292 


13155(4) 


s 


93497 


52872 


13709(5) 


c 


f2838 


52293 


13636(3) 


s 


93779 


52874 


16180(6) 


c 


f6596 


52320 


11289(4) 


s 


97951 


52904 


10985(5) 


s 


f683f 


52321 


79975(7) 


s 


99892 


52918 


04294(5) 


s 


f6968 


52322 


78333(5) 


s 


134317 


53165 


15846(4) 


s 


f7107 


52323 


78128(5) 


s 


134875 


53169 


16426(7) 


s 


43614 


52514 


05836(2) 


s 


146279 


53251 


02627(6) 


s 


43745 


52514 


99882(4) 


s 


146847 


53255 


10348(3) 


s 


44162 


52517 


99198(3) 


s 


167703 


53404 


81578(7) 


c 


47844 


52544 


42274(3) 


s 


167841 


53405 


80629(6) 


c 


47982 


52545 


41327(2) 


s 


167981 


53406 


81140(8) 


c 


48121 


52546 


41115(2) 


u 


196318 


53610 


22484(1) 


c 


48392 


52548 


35663(4) 


u 


196454 
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20136(1) 


c 


48675 
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38/96(2] 


u 


196595 
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zl349(lj 


s 


48763 
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u 


197006 


53615 


16357(2) 
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49035 
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97227(4) 
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245612 


53964 


07579(4) 


s 


49051 


52553 


08725(4) 


s 


245762 


53965 


15266(7) 


c 


49096 


52553 


41030(3) 


c 


259428 


54063 


25214(5) 


s 


49158 


52553 


85520(5) 


c 


265117 


54104 


09019(8) 


s 


49181 


52554 


02031(3) 


c 


315130 


54463 


10209(8) 


c 


49235 


52554 


40809(4) 


c 


317361 


54479 


11720(6) 


s 


49302 


52554 


88571(38) 


c 


319869 


54497 


12051(5) 


s 


49441 


52555 


88667(4) 


c 


359186 


54779 


35267(4) 


c 


49509 


52556 


37485(2) 


s 


405672 


55113 


04683(4) 


s 


52360 


52576 


84022(5) 


c 


417121 


55195 


23200(8) 


s 


52535 


52578 


09646(8) 


c 


417387 


55197 


14168(5) 


s 


54316 


52590 


88107(6) 


u 


418324 


55203 


86763(2) 


s 


56127 


52603 


88493(28) 


A 


449996 


55431 


22123(2) 


s 


56682 


52607 


86516(5) 











unclear, but it is most likely associated with the accretion 
disc. 



4 ORBITAL TIMINGS 

In order to measure the photometric period in ES Get we 
choose to fit a simple function to each night of data. We 
found that the variation in ES Get is fit well on any given 
night by the sum of three sinusoids of the form f{t) — 
a„ cos(27™(f — to — tn)/P) where n is 1, 2 or 3, hence these 
sinusoids correspond to the fundamental orbital frequency 
and the first and second harmonics. P is the orbital period 
and to is the time of minimum light. We fitted this function 
to each night of ULTRAGAM g'-hand data (in which the 
signal-to- noise is highest) . These fits are over-plotted on the 
g'-band data in Figure [1] We first set to to and P to 1 
and fit the function to the phase folded light curve in or- 
der to determine a„ and t„. We then determined to from 
the unfolded light curve by fixing a„ and t„ and fitting the 
function to these data. The process of fitting the data is 



complicated by the variable shape of the light curve, which 
we discussed in Section [31 This particularly complicates the 
task of fitting the GBA/SAAO observations, which have a 
much lower cadence than our WHT observations. It is there- 
fore difficult to fit most of the observations with the same 
degree of precision with which we fit the WHT observations. 

Furthermore, we find that our determination of the light 
curve minimum to is perturbed by the changing light curve 
shape. This is apparent in the 2005 ULTRAGAM timings. 
These were collected over the course of a few successive 
nights, and so when the linear ephemeris is subtracted they 
should have similar O ~~ C values. We find the O — C val- 
ues to form two groups with a significant offset of around 
20s (with a formal error on the individual timings of ~ls). 
These two groups correspond to the two distinct light curve 
shapes observed in Figure [1] In an attempt to characterise 
this effect with the aim of subtracting it from our light 
curves, we fitted each night separately and examined the 
variation in the model parameters with the O — C excess 
from a linear ephemeris. We found first of all that the model 
parameters can take a continuous range of values between 
the two extreme states we observed with ULTRAGAM, and 
we observed a scatter of up to ~25s in the timings. Some 
component of this scatter is doubtless due to flickering, the 
stochastic variation observed in all accreting systems; but 
a large component is due to the gross changes in the light 
curve. We observed a loose correlation between the O — C 
values and the amplitude and phase of the first harmonic. 
We attempted to improve our O — C residuals by adding 
an offset determined from these correlations, but since the 
correlation is fairly weak our timings were not improved to 
any significant degree. The limitation here is most likely the 
difficulty in fitting the majority of our observations with 
precision, due to the low time resolution and the gross vari- 
ation of the light curve on timescales of a few cycles. Since 
any improvement is very small, we chose not to apply any 
correction of this nature to our final determinations of the 
orbital timings. 

For consistency we obtain our final timings by applying 
the same model fit to each night, allowing to to vary but 
fixing an and tn- In practise we found it does not matter 
which night of ULTRAGAM data we use for these a„ and 
tn values, but we choose night 6 since the light curve shape 
in this case is intermediate between the two extremes. We 
did try phase-folding and fitting each night separately, and 
then using the fit for each night to determine the time of 
minimum light for that individual night, but this did not 
improve the timings to any significant degree and in some 
cases made things worse due to the uncertainty of the fits. 
We obtain 65 timings in total, which we list in Table [T] Our 
final set of timings is in good agreement (RMS ~6s) with the 
measurements of Espaillat ct al. (2005) for the cases where 
we share data. We convert all timings to the barycentric 
dynamical timescale, correcting for light travel times. 

Finally, using these timings, we calculate an ephemeris 

of 

BMJD{TDB) = 52200.980575(6) + 0.00717837598(3)S 

for the minimum-light point of the orbital cycle, using 
a linear least-squares fit. The data points in Figure [2] 
show the residuals after this linear fit is subtracted, and 
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demonstrates that there is no significant departure from 
linearity over the baseline of our observations. 



5 DISCUSSION 

In this section we examine the period change in ES Get over 
the baseline of our observations. In Section 15.11 we discuss 
the period changes and subsequent departure from a linear 
ephemeris in AM CVn systems as a result of gravitational 
radiation. We then compare these models to our ES Get 
timings in Section [5.21 



5.1 Period changes in accreting binaries 

The fractional period change in an accreting binary system 
can be shown to be 



(1) 



iMarsh fc Nelemand 120051 ) assuming conservative mass 
transfer, and where q = M2/M1, Mi and M2 are the masses 
of the accretor and donor respectively, and M2 is the mass 
transfer rate. The corresponding quadratic time shift in the 
times of light curve minima is given by 



AT = 



^2 



(2) 



where ts is the baseline of our observations. Assuming that 
the spin of the acc retor is strongly coupled to the orbit (see 
iMarsh et al.ll2003 ). the J term in Equation [T] is simply due 
to gravitational radiation, where 



J orb 



32 G-^ M1M2M 



(3) 



IILandau fc LifshitdHoTSI ). where a is the orbital separation. 
Since J/ J < in the absence of mass transfer, we would 
expect a decreasing period as the two components spiral 
in. In an accreting system the M2 term in Equation [1] can 
compensate for J (since M2 < 0) and for the expected mass 
transfer rates in AM GVn stars it is larger, and so we expect 
a lengthening period in these systems once it has settled into 
its long-term, stable state. It can be shown that 



J/J 



M2_ 

M2 ^ 1 + (C2 - Crl)/2 - q 
l|Marsh et al.ll2004l ). where 
d log R2 



C2 = 



and 



Crl 



d log M2 ' 



dlog(i^L/a) 
dlogM2 



(4) 



(5) 



(6) 



describe the response of the donor and the Roche lobe to 
mass transfer. Crl is taken to be 1/3 using the small q ap- 
proximation of IPaczvhskH (|l97ll ) . Since the response of a 
degenerate star is to expand on mass loss, ^2 is < 0. A value 
of = —0.28 is appropriate for a fully degenerate donor, 
but as the degre e of degeneracy of the donor is decreased, 
^2 approaches (|Delove et al.ll2005l ). 



5.2 The period evolution in ES Get 

In this section we compare our measurements of the orbital 
timings with the predicted time shift over our baseline. We 
see in Section [5.11 that the predicted time shift depends on 
the masses of the donor star and the accretor. We begin 
in Section 15.2.11 by selecting various appropriate values for 
these unknown parameters. We then compare our timings 
with the expectations for these parameters in Sections 15.2.21 
and 15.2.51 We consider first the case in which there is mass 
transfer between the two components. We also consider the 
detached case, in which there is no mass transfer and the pe- 
riod decreases at the rate given by Equation [T] with M2 — 0. 
However given that the spectral observations imply accre- 
tion, we would expect the first case to be a better fit to our 
data. 



5.2.1 The component masses 

There are three formation paths proposed for AM GVns. 
All three paths are consistent with a donor that is par- 
tially degenerate to some degree, and measurements of the 
component masses in other AM CVn syste ms have con- 
firmed the donors to be partially degen erate (|Roelofs et al.l 
l2006l . I2OO7I: iGopperwheat et al.l 2010 bl) . The 'white dwarf 
channeT (|Nelemans et al.l I2OOII ) suggests detached close 
double white dwarfs which are brought into contact as 
a result of angular momen tum loss due to gravitational 
wave radiation fGWR). Nelemans et al.l (120011 ) used a fuUy- 
deg enerate, zero- t empe rature donor in their formulation, 
but iDelove et al.l (I2OO5I ) argued that the donors could be 
semi-degenerate depending on the contact time of the bi- 
nary. The second scenari o is the 'evolved GV channel' 
l|Podsiadlowski et al.ll2003l ) which suggests the progenitors 
of AM GVns are GVs with evolved secondaries. The donor 
in this channel is initially non-degenerate and hydrogen-rich, 
but becomes degenerate and helium-rich (but still with a 
few per cent hydrogen) during its evolution before Roche 
lobe overfiow. The third sc enario is the 'helium star chan- 
nel' l|lben fc Tutukovlll99ll ). in which the donor had been a 
helium core burning star before coming into contact with its 
Roche lobe. For the white dwarf and evolved GV channels 
the donor mass is somewhat arbitrary, since it depends on 
factors s uch as the specific entropy of the donor at time of 
contact (|Delove et al.ll2005l ). We therefore chose to examine 
two cases: the first in which the donor is fully degenerate, 
and the second in which it has a mass consistent with a he- 
lium star progenitor. The helium star channel implies a very 
massive donor star, and the expected mass for the other two 
channels will lie between these two extremes. 

If we ass um e a fu lly degenerate donor, we find A'h to 
be 0.062Mf7, (iDelove et al.ii2005l ). This is the lowest possible 
mass for the donor in this system. We chose to combine this 
with the lowest likely accretor mass Mi, since the case in 
which both Mi and M2 are at their lowest values produces 
the minimum period change in both the accreting and non- 
accreting formulations. For this we assumed the system is 
accreting via a disc, rather than through the direct impact of 
the accretion stream onto the white dwarf primary. We be- 
lieve accretion via a disc is a reasonable assumption, since 
the spectroscopic data do feature the double peaked lines 
which are characteristic of an accretion disc, although the 
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direct-impact scenario cannot be ruled out entirely, since 
even in this case the primary Roche lobe might be filled 
with diffuse, hot gas which could produc e a 'disc-like' sig- 
nature in the spectra (jRoelofs et al.ll2010l ). The presence of 
an accretion disc implies a minimum mass for the accretor 
(for a given donor mass) simply via geometric arguments, 
since the accretion stream follows a ballistic trajectory. For 
a donor mass of O.O62M0, the minimum possible mass for 
the accretor is QAAMq. This accretor mass is quite low, so 
we also considered a higher accretor mass of ^Q.QMq, which 
is more typical of field white dwarfs. Since the donor is fully 
degenerate, we set 1^2 = —0.28. 

For the partiall y degenerate case, we use the 
Nelemans et al.l (l200ll }_ fit to the He-star M-R track of 



Tutukov fc Fedoroval ( 1989) to determine the donor mass. 
Figure 1 of iDelove et a l. (2005) shows that the predicted 
donor mass for this scenario is ~0.26Mq. We again assume 
the accretor mass to be at the minimum for accretion via a 
disc, which for this donor mass is 0.69Mq . The value of ^2 in 
this case is not clear but will lie somewhere between and 
the fully-degenerate value of —0.28. Given that the donor 
star we consider here is considerably less degenerate than 
the M2 = O.O62M0 case we would expect ^2 to be closer to 
0, but we consider both of these two extremes. 

5.2.2 Period evolution with mass transfer 

We plot in Figure [2] the deviation from a linear ephemeris 
for our orbital timings in terms of phase, using the linear 
ephemeris given in Sectional If we examine the data points 
first of all, we see a large scatter (as discussed in Section |4]) 
is apparent in these timings (in particular, note the separa- 
tion between the two groups of 2005 ULTRACAM timings, 
at ~200, 000 cycles). However, aside from this random vari- 
ation, there is no obvious departure from a linear ephemeris 
over the baseline of our observations. There might be a very 
small increase in period over our baseline (manifested by a 
slightly positive O—C value at either extreme of the plot and 
a slightly negative value in the middle), however any such 
trend does not appear to be very significant when compared 
to the noise. 

The curved lines in the top panel of Figure [2] show 
the predicted departure from a linear ephemeris with 
Ml = O.44M0, M2 = O.O62M0 (solid line); Mi = O.6M0, 
M2 = O.O62M0 (dashed hue) and Mi = O.69M0, M2 = 
O.26M0(the dot-dashed and dotted lines, for a ("2 of —0.28 
and respectively). All four predict an increasing period 
as a result of the mass transfer. These lines have been ver- 
tically offset to optimise the fit to the data. We see that 
in the cases where we assume a fully degenerate donor the 
predicted time shift is very small, and the amplitude is less 
than the uncertainty in our data points. A fully degenerate 
donor is therefore consistent with our timings. Increasing 
the mass of the accretor has little effect on these ampli- 
tudes. Conversely, the predicted time shift when a donor 
mass of A/2 = O.26M0 is used is much larger. When we set 
C,2 to —0.28 it is very large, and clearly precluded by our 
measurements. As ("2 is increased the response of the donor 
to mass loss is decreased, and so the expected time shift 
over our baseline is also decreased. The exact value of (,2 for 
this semi-degenerate donor is likely close to (but not greater 
than) 0. The dotted line in Figure [2] shows that even when 



we set C,2 = Q the predicted deviation from linearity is not 
consistent with our timings if we assume a massive donor. 

We have shown that the He-star formation channel pro- 
duces a donor that is too massive to be consistent with 
ou r timings. If we examine t he evolutionary tracks given 
in IPodsiadlowski et al.l l|2003l ) for the evolved CV channel 
we find that, while the donor mass is dependent on initial 
conditions, a reasonable assumption for the donor mass is 
0.11 - O.I2M0 , which also predicts a deviation from linearity 
which is inconsistent with our timings. We conclude there- 
fore that the double white dwarf channel is the most likely 
formation path for the low donor mass we infer. 

An alternative explanation is that ES Get is in the ini- 
tial, turn-on phase before the long AM GVn phase of length- 
ening period, and so the accretion rate in ES Get is currently 
below the stable rate which would be expected if it were in 
that phase. The decreasing periods of HM Cue and V407 
Vul have been explained by invoking an accretion model in 
which the donor has a thin hydrogen envelope, and it is 
the acc retion of this env elope which leads to the shortening 
period (|D'Antona et al.|[200S V At some point this envelope 
will be depleted and the transfer of helium will begin, leading 
to a lengthening period along the 'traditional' evolutionary 
track. As the binary transitions from one state to the other 
the fractional period change must drop significantly. Since 
we find ES Get to be consistent with a P of it could now 
exist in this intermediate regime, linking HM Cue and V407 
Vul with the longer period AM CVns. 

Finally, we note that one assumption we have made is 
that all of the angular momentum contained in the mass 
gained by the accretor from the donor is fed back into the 
orbit. Alternatively, some fraction of the disc accretion may 
act to spin up the accretor. However, we find that this only 
increases the expected P, and so does not improve the fit 
of any of our models to the observed timings. If we consider 
the extreme case, in which none of the angular momentum 
is fed back to the orbit. Equation |4] is modified to 



M2 
M2 



J/J 



1 + (C2 - CiiL)/2 - g - + 



(7) 



(jMarsh et al.ll20o3 ). where Ri is the accretor radius and a 
is the binary separation. This extreme case predicts a devia- 
tion from linearity of 0.23 or 0.32 in phase over our baseline, 
assuming a fully degenerate donor and an accretor mass of 
0.60 or O.45M0. This is clearly inconsistent with our tim- 
ings, suggesting that a substantial fraction of the angular 
momentum is indeed fed back into the orbit. 



5.2.3 Period evolution with no mass transfer 

We consider also the scenario in which the donor star is de- 
tached from its Roche lobe. In this case the period will de- 
crease as the components spiral in as angular momentum is 
lost due to gravitational radiation. The three curved tracks 
in the bottom panel of Figure[2]show the predicted departure 
from a linear ephemeris for the detached case, for the same 
combinations of Mi and M2 used in the top panel. All three 
lines are a poorer fit to our timings compared to the cor- 
responding semi-detached scenarios, as expected given the 
evidence of accretion. 
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Figure 2. Top panel: Minimum-light timings for our complete dataset. We plot the period excess (O — C) in terms of phase against 
the orbital cycle. ULTRACAM observations are denoted by the black crosses, SAAO observations by the unfilled red circles, CBA 
observations by the unfilled blue squares, and the ACAM observation by the filled grey square. The four curved lines indicate the 
expected period change when we include mass transfer at the secular rate, for different combinations of accretor and donor mass. The 
four combinations are (1) solid line, M\ = OAAMq, M2 = 0.062Mq, (2 = —0.28. This is the minimum period change, assuming accretion 
via a disc and a fully degenerate donor, (ii) dashed line. Mi = O.6OM0, M2 = 0.062Mq, 1^2 = —0.28. This is a variation on the minimum 
case in which we assume a more likely accretor mass, (iii) dot-dashed line, Mi = O.GQM q, M2 = 0.26M(7), C2 = —0.28. Here we assume 
a partially degenerate donor using the M-R relation for the He-star formation channel l lNelemans et aUl200ll ). and the corresponding 
minimum accretor mass. We use the same value of ^2 as for the fully degenerate models. This is unlikely to be correct, and so in (iv, 
dotted line) we set ^2 = 0, keeping all other parameters the same. The true value of (^2 will be between these two extremes, but is likely 
to be closest to 0. Bottom panel: As with the top panel, except here we assume the secondary star is detached from its Roche lobe, 
and so the binary evolution proceeds due to gravitational radiation in the absence of mass transfer. In this detached case models (iii) 
and (iv) give the same prediction, since they are only distinguished from each other by the response of the donor to mass loss. 



6 CONCLUSIONS 

We have obtained high time resolution lightcurves of the 
AM CVn star ES Get with WHT-hULTRACAM and other 
smaller telescopes around the world. We observed the previ- 
ously reported cyclical variability which has been confirmed 
to be orbital in origin. Our observations cover a baseline of 
9 years and we fit the lightcurves to obtain precise orbital 
timings over that period. We find that there is no evidence 
for a deviation from a linear ephemeris over this time period, 
although scatter in the timings is high due to fiickering and 
large cycle-to-cycle variations in the hght curve shape. There 
may be an underlying trend which is masked by this scatter. 
However, many predictions do suggest a lengthening orbital 
period which should be clearly detectable over this noise. 
Our findings suggest a low mass for the donor star, close or 
equal to the zero temperature mass of 0.062AfQ. This would 
make the double white dwarf formation channel the most 
likely scenario for this system. Alternatively, the accretion 
rate in this system may be significantly below the long-term. 



secular rate for an accreting AM CVn star. ES Get is the 
shortest period AM CVn star and the only ultracompact bi- 
naries with shorter periods are HM Gnc and V407 Vul, both 
of which have been shown to have a decreasing period. ES 
Get may be in an intermediate phase between these systems 
and the longer period AM GVns. 
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